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We test an algebraic algorithm based on the coordinate-space method, evaluating with high accuracy the critical 
mass for Wilson fermions in lattice QCD at two loops. We test the results by using different types of infrared 
regularization. 



We have already presented at a Lattice confer- 
ence |j| an algebraic algorithm that allows to ap- 
ply the coordinate-space method by Liischer and 
Weisz {I to two-loop lattice integrals with gluon 
and Wilson-fermion propagators. In order to test 
the method, we have recently redone H a two- 
loop computation of the critical mass for Wilson 
fermions 0. 

For Wilson fermions, the dressed inverse 
fermion propagator has the form 

S~ 1 (p,m B ) = 

= if +m B + M w (p) - £ L (p, ms^go), (1) 

where, setting the lattice spacing equal to one, 

= sinp^ , (2) 

f = J2 (2 sin |) 2 , (3) 

Mw(p) = r -yP 2 ■ (4) 

The additive mass renormalization Stub is ob- 
tained by requiring S'~ 1 (0, Stub) — 0, i.e. 

S L (0, Sm B ,g ) = Sm B - (5) 

This equation can be solved in perturbation the- 



ory by expanding 

oo 

£ L (0,m B , 5o H]T 5o 2 "£("). (6) 

71=1 

We have computed £W and S^ 2 ^ for = 1, 
gauge group SU(N), and Nf fermionic flavour 
species, in the Feynman gauge. 

In Ref. H we already reported the analytic one- 
loop expression for the fermionic self-energy E L . 
Our result was expressed in terms of three purely 
bosonic constants Zq, Z\ and Fq and of 12 nu- 
merical constants that appear in the presence of 
Wilson fermions. The numerical values of these 
constant are obtained by using a powerful recur- 
sive method that gives very precise results ||,[|. 
This algorithm generalizes the method we intro- 
duced for purely bosonic integrals in j?J . In prac- 
tice, we have computed all constants but Fq with 
60-digit precision. 

At one-loop order 

i=l 

where are the contributions of the two con- 
tributing diagrams. In terms of the basic integrals 



2 



they are given by 



= -Z 

= f-^(i,o) 



Summing up the two contributions we obtain 
2 



E- 



[Z + 2^(1,0)] 
0.16285705871085078618 . (8) 



The constant is in excellent agreement with the 

result of Ref. (§, 

^ 4 cf } = -0.162857058711(2). 

At two loops there are 26 diagrams. They are 
numbered as in Ref. Q in order to simplify the 
comparison. The i-th diagram gives a contribu- 
tion of the form 



Di = (N 2 - 1) 



i 2 ) 



1 

JV2 



..(2) 



Nf_ (2) 

N °^ 



(9) 



In Table | we report the results given in 
Ref. Q and those obtained here by using the 
configuration-space method. When we have not 
reported an error, the precision we achieve is 
much higher than the reported digits. This oc- 
curs in general when the diagram is the product 
of one-loop integrals. All results are in agree- 
ment with those presented in Ref. ||. Only for 
diagram 6 there is apparently a (very) small un- 
derestimation of the error, which is negligible in 
the sum of all contributions. In Table ^ diagrams 
are grouped together in order to obtain infrared- 
convergent results, but this not necessary in our 
method. Indeed, we can compute each of them 
separately, by introducing an infrared regulariza- 
tion. To test the results, we have used four dif- 
ferent infrared regularizations: 

(a) We add a mass in the denominators of 
the propagators. Explicitly, for the gluon 
(As(fc)) and for the fermion (Ap(k)) prop- 
agator we use: 



1 



fc 2 



t + M w (k) 2 +m 2 



(b) We regularize the gluon propagator as in 
(a), but use instead the correct Wilson- 
fermion propagator 



A F (fc) = 



-ik^ + M w (k) + m 
k 2 + (M w (k)+m) 2 



(c) We regularize the Wilson fermion as in (a), 
but use the massless propagator for the 
gluon. 

(d) We regularize the Wilson fermion as in (b) 
and the massless propagator for the gluon. 

As an example, we report in Table || the diver- 
gent and the finite contribution of three diagrams 
whose sum is infrared finite. If we write 



D 



(x) 



D 



(a) 



1 

iV 2 " 



A 



(x) 



where (x) refers to the chosen infrared regulariza- 
tion, we get the results reported in Table ||. Indi- 
vidual diagrams depend on the regularization but 
their sum does not. 
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c< 2 > 






i,i 


2,i 


3,i 
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0.002000362950707492 


-0.00300O.-, 1 1 1200612375 







0.0020003629507074987148 


-0.0030005444260612480722 





4 


0.00040921361(1) 


-0.00061382041(2) 







0.0004092136068803147865 


-0.0006138204103204721798 
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6 


-0.0000488891(8) 


0.000097778(2) 







-0.0000488853119(2) 


0.000097770623(5) 





7+8 + 9+10 + 11 


-0.013927(3) 


0.014525(2) 







— 0.01392647740 (2) 


0.0145250053341618950704 


o 


12 + 13 








0.00079263(8) 










0.000792647(2) 


14+15 + 16+17+18 


-0.005753(1) 


0.0058323(7) 







-0.00575248584(1) 


0.005832127004694453 





19+20 








0.000393556(7) 










0.000393556(4) 


21 + 22 + 23 


0.000096768(4) 


-0.000096768(4) 







0.0000967648(2) 


-0.0000967648(2) 
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25 


0.00007762(1) 


-0.00015524(3) 







0.000077613106(4) 


-0.000155226212(8) 





26 


-0.00040000(5) 










-0.00039997586(1) 








27 





-0.000006522(1) 










-0.0000065203(1) 





28 


0.0000078482(5) 


-0.000015696(1) 







0.000007S iMM.r, 2 72 2033294 


-0.0000156961305444066589 





Total 


-0.017537(3) 


0.016567(2) 


0.00118618(8) 




-0.0175360218(2) 


0.0165663304(2) 


0.001186203(6) 



Table 1 

Coefficients c{\ , Sf\ and . For each of them we report in the first line the result of Ref . Q , obtained 
by means of a momentum-space integration, and in the second line our result, obtained by means of the 
coordinate-space method. 



~i Divergent Part c^] 
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(1- 






9 

;m 




0.001606284825541242 


-0.001606284825541242 
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(1- 






JF(1,0) > 


| log TO 2 


0.0005015205(2) 


-0.0005015205(2) 





23 


(1- 






^(1,0)' 
8tt 2 , 


| logm 2 


-0.002011040454066014 


0.002011040454066014 






Table 2 

Divergent and finite contribution of three diagrams whose sum is infrared finite. 
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T 96tT 2 




~ 96tT 2 
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32tT 2 


^(1,0) 

16tT 2 


llZn 
192ir 2 


, 11^(1,0) 

96tT 2 


16Z 

1927T 2 


. 16^(1,0) 
+ 96tt 2 
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z 

32?r 2 


^(1,0) 

16tT 2 


llZo 

1927T 2 


11^(1,0) 

967T 2 


16Z Q 

1927T 2 


16^(1,0) 
96tt 2 



Table 3 

Contribution from different regularizations. 



